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Semi-convection: What is the underlying 
physical context ? 


A. Noels 


Abstract Stellar conditions leading to a possible semi-convective mixing are dis¬ 
cussed in three relevant cases: (1) low mass MS stars in which the CNO cycle takes 
progressively the lead over the PP chain due to the increase in temperature as core 
hydrogen burning proceeds, (2) massive MS stars which experience a large contri¬ 
bution of the radiation pressure to the total pressure and (3) core helium burning 
stars for which the production of carbon in the core increases the opacity. A short 
discussion of semi-convection in terms of instability of non radial modes follows. 


1 Introduction 

In main sequence stars massive enough to burn hydrogen through the CNO cycle, a 
convective core is already present at the ZAMS. In most cases, its mass extension is 
maximum at the ZAMS and then it shrinks with time due to the decrease in opacity 
resulting from the transformation of hydrogen into helium which rather drastically 
reduces the density of free electrons. In such a case, the opacity is larger outside 
the convective core in layers which are richer in hydrogen but as a result of the pro¬ 
gressive decrease in mass extension of the convective core; there is neither chemical 
nor opacity discontinuity at the convective border. In some cases however the con¬ 
vective core mass tends to grow as the star evolves and this leads to the formation 
of a hydrogen discontinuity at the border, the outer border being hydrogen richer 
and thus more opaque than the inner border. The layers affected by this, although 
outside the convective core, are unstable towards convection if the Schwarzschild 
criterion 03 is used to define the convective neutrality, i.e. 
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where the ^rad and refer to radiative and adiabatic temperature gradients and Pi 
is the secong adiabatic coefficient, while they are stable when applying the Ledoux 
criterion a since 
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where )3 is the ratio of the gas pressure to the total pressure and /r is the mean 
molecular weight. In a pioneering analysis, Schwarzschild and Harm M solved 
this problem by adding a partial chemical mixing in the so-called semi-convective 
layers in order to ensure their convective neutrality. The helium abundance in these 
layers was enriched to just the precise amount required to satisfy the Schwarzschild 
criterion. A similar partial mixing was adopted by Sakashita and Hayashi lITSll but 
with the Ledoux criterion instead. 

After recalling a few basic points of stellar structure in section]^ I shall discuss 
the physical conditions leading to such an increase with time of the convective core 
mass, first in low mass stars (sections]^ and then in massive stars (section]^. In 
section]^ I shall briefly address the problem of semi-convection in helium burning 
stars. Section |6] will be devoted to a short discussion in terms of vibrational sta¬ 
bility. My aim is here to emphasize the physical conditions for semi-convection to 
(possibly) develop, not to describe the modem ways of tackling the problem nor to 
present an exhaustive review of the theoretical works done since the sixties. These 
are presented and discussed in the next chapter US of which this constitutes a sort 
of preamble. 


2 A few basic points 

Let us recall here a few basic points affecting the stellar structure ; 

1. Radiative temperature gradient 

In a simplified way the radiative temperature gradient can be written 

'^rad^-K (3) 

m 

where L is the luminosity and K the opacity. This means that 

• a large L/m value typical of nuclear burning cores is favorable to convection. 
The larger the temperature sensitivity of the nuclear energy production rate, 
the larger the L/m value; 

• a large opacity mostly found in ionization zones located in the outer layers 
leads to the presence of a convective envelope. 

2. Temperature sensitivity of a nuclear reaction 
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The temperature sensitivity of a non resonant nuclear reaction involving the fu¬ 
sion of two nucleons and (see for instance 0) is given by 


v = (^) 


( 4 ) 


where e is the nuclear energy production rate, T is the temperature and p the 
density. Its value strongly depends on the Gamow factor b 


b = ZaZ„A 


1/2 


with 
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The V value is then easily estimated from the relations 


V = 


t-2 


with 


T = 


3 ^ 

kT 


and 


= (6) 


where k is the Boltzmann constant and Eq is the Gamow energy Le. the most 
effective energy for the nuclear reaction to take place. 


3. Temperature sensitivity of the PP chain reactions 

From the relations above it is evident that the smallest v value will be obtained 
for the reaction. At a temperature of about 10^K (kT ^ 1 keV) Vn is of 

the order of 4. As soon as the charge of the nucleons increases as in (^He,\He), 
the Gamow factor b drastically increases and so does the temperature sensitivity. 
For kT equal to 1 keV, V33 reaches a value of about 17 . 

• PP chain operating out of equilibrium 

Near the end of the pre main sequence phase when the temperature at the cen¬ 
ter reaches a value of about 10^K the hydrogen burning nuclear reactions start. 
In low mass stars hydrogen burning is largely dominated by the PP chain. The 
abundances of the nucleons involved in the PP chain, essentially that of 
are however still different from their equilibrium (or more precisely station¬ 
ary) abundances. This means that \He is accumulating up to the point where 
its destruction rate will equal its formation rate. The temperature sensitivity is 
given by 

£11 , £33 

Vpp = Vi 1-h V33 — ( 7 ) 

e e 

Due to the high value of V33 the resulting sensitivity is large and a convective 
core appears and remains during the whole process of reaching the equilib¬ 
rium abundances. 

• PP chain operating at equilibrium 

When \H and reach their equilibrium abundances i.e. when they are pro¬ 
duced and destroyed at exactly the same rate, the whole PP chain is governed 
by the ([//,[//) reaction (or pp reaction). That means that the v value is low 
and the L/m ratio is not high enough to allow the presence of a convective 


core. 
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• PP chain operating with overshooting 

If a certain amount of overshooting, or any other extra mixing, is taken into 
account above the convective core boundary, the convective core present be¬ 
fore ZAMS can be maintained during the whole main sequence phase. This 
is due to the fact that this extra mixing prevents \He from reaching its equi¬ 
librium abundance whatever the elapsed time since fresh \He is continuously 
brought into the core. 

4. Temperature sensitivity of the CNO cycle reactions 

Whatever the proton capture reaction involved in the CNO cycle, its sensitivity 
is high due to the large value of the Gamow factor. A convective core is present 
during the whole core hydrogen burning phase. 

5. Temperature sensitivity of the 3 a helium burning reaction 

The 3 a reaction is a resonant reaction. Its temperature sensitivity (see 0) is 
given by 

42.9 

V3a = ^-3 (8) 

where Tg is the temperature in lO^/T. Such large values of v involve the presence 
of a convective helium burning core in the whole stellar mass domain. 

6. Central temperature of MS stars 

Assuming hydrostatic equilibrium, thermal equilibrium and radiative transfer, a 
dimensional reasoning easily leads to (see for instance 0) 

M 1 v-i 

T -(9) 

R 

Where M is the stellar mass and R is the surface radius. Whatever the way of 
burning hydrogen, through PP chain or CNO cycle, the exponent of M is positive 
which means that the central temperature increases with the stellar mass. This 
implies a progressive growth of the contribution of CNO in the nuclear reactions, 
from PP chain for low mass stars to CNO cycle in intermediate and massive 
stars. Another important effect is the increase of the contribution of the radiation 
pressure to the total pressure as the stellar mass increases. 


3 Semi-convection in low mass stars 

In the mass range [I.QMq — 2.QMq] the small convective core present while the PP 
chain elements are still reaching their equilibrium values, vanishes on the ZAMS. 
As the star evolves however the temperature at the center slowly increases which 
leads to a growing contribution of the CNO cycle to the nuclear reactions. The tem¬ 
perature sensitivity increases from a value of about 5 to a larger value typical of 
CNO reactions (~ 15). A convective core appears and its mass extension increases 
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as the CNO contribution increases. Figure shows the evolution of the fractional 
core mass extension with the central hydrogen abundance. This growing tendency 
is enhanced when models are computed with overshooting since a convective core 
is already present at the ZAMS due to the non equilibrium values of the \He abun¬ 
dance in the mixed region (see section]^. 



Fig. 1 Fractional convective core mass as a function of central hydrogen abundance Xc for models 
in the mass range 1 .OMq 2.OM0 computed without and with overshooting. The overshooting pa¬ 
rameter jS (ratio of the overshooting distance and the minimum between the convective core radius 
and the local pressure scale height) is indicated in each panel. (From GHl) 


When models are computed with the Ledoux criterion such an evolution of the 
convective core leads to a well dehned and numerically stable discontinuity at the 
convective core boundary as can be seen in the left panel in Fig. This is not 
the case when the Schwarzschild criterion is used (see the right panel in Figure 
1^ since small convective shells form in the region of varying mean molecular 
weight (/f-gradient region). The resulting effect is to create nearly neutral convec¬ 
tive conditions in the whole /f-gradient region and to reproduce more or less the 
Schwarzschild and Harm solution of the problem of semi-convection (see section 
[T]i. The number and the extent of these small convective shells are however depen¬ 
dent on the number of mesh points in the model in a somewhat erratic way. Although 
rather tempting this solution is dangerous especially because of the potential numer¬ 
ical diffusion induced by these moving convective layers. The mass domain of stars 
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affected by this semi-convective problem is limited to [1 .OMq — 2.0Mq] and the du¬ 
ration of the phase itself is only a small fraction of the main sequence lifetime since 
those intermediate convective zones appears mostly near the maximum extent of the 
convective core. 



Fig. 2 Hydrogen profile in models of I.SMq (X = 0.70, Z = 0.02) computed with the Ledoux cri¬ 
terion (left panel) and with the Schwarzschild criterion (right panel) in the fractional mass interval 
[ 0 , 0 . 2 ]. 


For intermediate mass stars more massive than 2.0Mq the problem of semi¬ 
convection disappears since hydrogen burning is dominated by the CNO cycle al¬ 
ready at the ZAMS and as the evolution proceeds the temperature sensitivity does 
not change significantly. Since the hydrogen abundance X decreases in the convec¬ 
tive core, V^ad decreases accordingly (see Figure]^ and the /r-gradient region forms 
a smooth transition between the convective core and the homogeneous envelope. 
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Fig. 3 Radiative and adiabatic temperature gradient distributions (left panel) and hydrogen profile 
(right panel) in models of SMq (X=0.70, Z=0.02) during MS. (From fTTI ) 


4 Semi-convection in massive stars 

Another problem arises for stellar masses larger than about 15Mq. As recalled in 
section the radiation pressure becomes a larger and larger contributor to the to¬ 
tal pressure as the stellar mass increases. The resulting effect is a decrease of the 
adiabatic temperature gradient due to a progressive decrease of i 2 from 5/3 to 4/3 
as j5 (see section]^ varies from 1 to 0. This effect favors convection and thus a 
larger mass extension of the convective core as M increases. As hydrogen burning 
proceeds the core mass extension decreases but and V^d are extremely close 
to one another in the /f-gradient region. With the Ledoux criterion the steepness 
of the /i-gradient is such that the layers remain stable towards convection. On the 
contrary when adopting the Schwarzschild criterion small convective shells appear 
in /f-gradient region. This is illustrated in Figure where the X profile, V^ad and 
Vad are shown in the left and the right panels for a typical model (Xc=0.20) of a 
MS 15M0 computed with the Ledoux and with the Schwarzschild criterions respec¬ 
tively. 
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Fig. 4 Hydrogen profile (dashed line), Vrad (full line) and (dotted line) for an MS model with 
Xc—O.lO) computed with the Ledoux criterion (left panel) and with the Schwarzschild criterion 
(right panel) in a sequence of ISMq MS (X=0.70, Z=0.02). 


Near the end of core hydrogen burning a smell convective zone appears at the 
base of the homogeneous envelope in models computed with the Ledoux criterion. 
This does not affect the /i-gradient profile which remains perfectly smooth with well 
defined junctions with the convective core and the homogeneous envelope plateaus 
during the whole MS. In models computed with the Schwarzschild criterion on the 
contrary the somewhat erratic behavior of transient convective shells leads to a step 
behavior of the X profile whose precise form depends on the adopted number of 
mesh points. The near equality of and Vad in the /r-gradient region tends to 
mimic the treatment of semi-convection advocated by m. 

As the mass increases [M > SOM©] the problem becomes more and more crucial 
since the mass extension of the convective cores even grows with time during the 
MS. The resulting discontinuity in X at the convective core boundary leads to a 
situation somewhat similar to that discussed in section|3] 


5 Semi-convection in helium burning stars 

In models of core helium burning massive stars the radiative opacity is dominated 
by electron scattering and is thus insensitive to the change in chemical composition 
produced by the transformation of helium into heavier a-elements. This is not the 
case however for models burning helium at lower temperature and higher density 
as it is found in intermediate and low mass stars. The free-free transitions opacity 
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depending on the charge of the ion is larger in carbon rich mixture than in helium 
richer ones. This is illustrated in the left panel of Figure|^taken from |[T|. 

Core helium burning stars all have a convective core as discussed in section 
As a result of the transformation of helium into carbon leading to a larger opacity 
the radiative temperature gradient increases in the convective core while it remains 
nearly constant just outside the convective boundary since the chemical composition 
is there unchanged. A discontinuity in tends to form at the core boundary as 
can be seen in Figure]^ where Vrad^ ^ad the helium profile as a function of the 

fractional mass [0,0.3] are given for a model of AModot (left panel) and SMq (right 
panel) in the core helium burning phase. 



Fig. 5 Helium profile (dashed line), (full line) and (dotted line) as a function of the 
fractional mass [0,0.3] in models of 4-Mq (left panel) and SMq (right panel) during core helium 
burning. A discontinuity in W^ad is visible at the boundary of the convective core. 


As discussed in m such a convective boundary is unstable since mixing a radia¬ 
tive layer close to the boundary with the homogeneous matter within the convective 
core enhances the carbon contain in the layer and makes it unstable towards con¬ 
vection. A sort of overshooting called by III a self driving mechanism leads to the 
extension of the core up to a layer for which the Schwarzschild criterion is satished 
with the chemical composition of the convective core. This progressive outward 
shift of the core boundary meets however a new difficulty which was first analyzed 
in 121. The mass distribution of indeed presents a minimum before reaching 
the boundary. Wether Vrad increases or decreases with time the existence of this 
minimum prevents a coherency in the determination of the convective boundary as 
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can be seen in the top and middle graphs in Figure]^ (right panel). The bottom graph 
shows the resulting induced semi-convection as proposed in la. 




Fig. 6 Ratio of opacities in a pure earbon and a pure helium matter as a function of the logarithm 
of the density for labeled values of the temperature expressed in 10*A" (left panel; from (T)). Mass 
distribution of and V„^ illustrating the problem of the minimum in Vrad during core helium 
burning. In the top and middle graphs V„,^ inereases and decreases with time respectively. The 
induced semi-convective mixing is shown in the bottom graph (right panel; from l2ll. 


6 Discussion in terms of vibrational stability 

I very briefly discuss here some earlier works on the stability aspects relevant to 
semi-convection. Classical convection is due to a dynamical instability of modes 
while overstable convection present in layers of varying mean molecular weight 
results from a vibrational instability of dynamically stable modes trapped in the 
/r-gradient region Q. The physical structures discussed in the above sections are 
generally favorable to the existence of such trapped modes. It is thus tempting to 
interpret a semi-convective mixing as a mixing resulting from the instability of these 
modes. The precise nature of the mixing is however far from straightforward. 
Since the destabilizing term tends to vanish when and V^d are equal, it would 
be safe to say that a mixing tending to this equality would lead to a stable situation. 
Another mechanism acting in massive stars was proposed by who showed that 
transient convective shells would progressively move across and partially mix the 
/r-gradient region until the Schwarzschild neutrality condition is met. 
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As already noted in section [T] attempts have however been made by numerous 
stellar evolution scientists at forming semi-mixed regions either neutral towards the 
Schwarzschild or towards the Ledoux criterion. As pointed by Ledoux El adopting 
the Schwarzschild criterion as a neutrality condition in a semi-convective region has 
the advantage of leading to a dynamically stable situation while the Ledoux criterion 
leads to marginally dynamically unstable conditions which would inevitably lead to 
a full convective mixing. Ledoux was indeed the first to advocate in this context 
the use of the Schwarzschild criterion and NOT the Ledoux criterion as a neutrality 
condition in semi-convective zones. 

Even trapped modes must be checked through a full stability analysis. In the case 
of massive MS stars Q showed that low order modes of high spherical harmonic 
degree could be trapped in the /i-gradient region. The timescales were found to be 
of the order of 10^ to 10"^ yr which made them good candidates for a partial semi- 
convective mixing. 

For low mass MS stars a similar analysis was performed and unstable trapped 
low order modes of high degree were also found. A very narrow mass range near 
1. IMq was however affected by this instability. Moreover unstable modes were only 
found during a short part of the main sequence phase near the maximum extent of 
the convective core. 

Although impossible to draw a full picture of semi-convection from these sta¬ 
bility analyses it is interesting to point out that they seem to converge towards a 
partial chemical mixing of the semi-convective layers ensuring their being neutral 
with respect to the Schwarzschild criterion. 
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